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ABSTRACT: High-pressure xenon gas is an attractive detection medium for a variety of applications 
in fundamental and applied physics. In this paper we study the transport properties of ionization 
electrons, and the mechanism of electron-ion recombination, in xenon gas at 10 bar pressure. For 
this purpose, we use a source of alpha particles in the NEXT-DEMO time projection chamber, 
the large scale prototype of the NEXT- 100 neutrinoless double beta decay experiment, in three 
different drift electric field configurations. Our electron drift velocity and longitudinal diffusion 
results are similar to expectations based on available electron scattering cross sections on pure 
xenon, favoring low-diffusion models. In addition, two types of measurements addressing the 
connection between the ionization and scintillation yields were performed. On the one hand we 
observe, for the first time in xenon gas, large event-by-event correlated fluctuations between the 
ionization and scintillation signals, similarly to what has already been observed in liquid xenon. 
On the other hand, we study the field dependence of the average scintillation yield. Both types of 
measurements may shed light on the mechanism of electron-ion recombination in xenon gas for 
highly-ionizing particles. 
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1. Introduction 

Particles traversing a xenon-based detector transfer energy to xenon atoms in the form of heat, 
ionization and excitation. The electron-ion pairs produced by the ionization process can give rise 
to a detectable signal, provided that a sufficiently strong electric field is applied to the detection 
medium to prevent electron-ion recombination. Scintillation photons emitted through the radiative 
decay of excited xenon molecules (or excitons) can also be detected. The possibility to detect both 
the scintillation and ionization signals produced by the passage of particles is the main advantage 
of xenon as a detection medium, together with its high stopping power. In addition, for xenon 
in the gas phase, the ionization fluctuations are characterized by a small intrinsic Fano factor flU], 
permitting accurate calorimetry using the ionization signal alone. For these reasons, xenon-based 
detectors are used for a variety of applications in fundamental and applied physics, including neu- 
trinoless double beta decay searches [0], WIMP dark matter direct detection QE ED, as well as X-ray 
astronomy, gamma-ray astronomy, and medical imaging 0. 0] . 
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The scintillation (often called SI)) and ionization (S2) signals measured by xenon-based de- 
tectors depend on the primary scintillation and ionization production yields, on the electron-ion 
recombination strength, and on the photon and ionization electron transport properties through the 
detector. The scintillation signal is also affected by the scintillation emission spectrum and time 
response. Because most xenon-based detectors make use of xenon in its liquid phase, the detection 
properties of liquid xenon [3D are generally better understood than the ones for xenon gas Q5D. A 
brief summary of xenon gas detection properties follows. 

The scintillation and ionization yields are usually expressed in terms of the average energies 
necessary to create a primary scintillation photon (W s ) and an electron-ion pair (WJ), respectively. 
Typically, both values are taken as empirical parameters provided by experiments. The mean en- 
ergy to form an electron-ion pair in high-pressure xenon, and when negligible recombination oc- 
curs, has been accurately measured to be about W\ = 22 eV using electrons and gamma rays [0, 0, |8D- 
The dependence of Wj on electric field and gas density has also been systematically studied. For 
gas densities lower than p ~ 0.3 g/cm 3 (gas pressures less than about 40 bars at room temperature, 
that is in the range of pressures of interest here), recombination effects are small and W[ does not 
depend on gas density [fj. For very high gas densities, p > 0.3 g/cm 3 , W{ depends not only on the 
electric field but also on the gas density fTTDTl . The understanding of the ionization yield in xenon 
gas for alpha particles is instead less complete. The much higher ionization density produced by 
alpha particles, together with more efficient recombination processes (see below), may quench the 
ionization signal. The ionization quenching factor q for alpha particles is an empirical value de- 
fined as the ratio of the amount of ionization produced by alpha particle to the amount of ionization 
produced by an electron (or gamma ray) for the same deposited energy. No ionization quenching 
(q = 1) is reported in [fH, where the value of W{ = 22.0 eV for electrons is to be compared with 
essentially the same number (Wi = 21.9 eV) for alphas. On the other hand, the authors of UTTH 
report a ionization quenching factor of q = 0. 15 in xenon gas at 5 bar and for a drift field of order 
1 kV/cm. 

The primary scintillation yield of xenon gas has traditionally been far less understood than the 
ionization yield. Recently, several measurements of W s have become available to try to remedy the 
situation. Value of W s = 72 ± 6 eV flI2, W s = 111 ± 16 eV flH and W s = 76 ± 12 eV [Q have 
been reported for xenon gas in the 1-20 bar range, using X-rays and low-energy gamma-rays. 

Electron-ion recombination affects not only the observable ionization signal, but also the scin- 
tillation one. As the recombination rate increases, the scintillation yield increases, given that each 
recombined pair is accompanied by the emission of one scintillation photon, sometimes called 
recombination luminescence. Recombination in high-pressure xenon gas has been studied by mea- 
suring how the average scintillation and ionization yields vary as a function of drift field and gas 
density, for both electrons and gamma rays [11311 as well as for alpha particles [1TH O] . For example, 
the authors of [|T7ll measure that the number of scintillation photons decreases monotonously with 
the electric field, and that the scintillation yield in the field-independent region is about 40% of that 
at zero electric field. In addition, in liquid xenon, recombination effects have also shown to pro- 
duce event-by-event, anti-correlated, fluctuations between the ionization and scintillation signals 
UTEl [TOT!. For example, the authors of HT9H report correlation coefficients between the ionization 
and scintillation signals in the —0.87 to —0.74 range for 662 keV gamma rays, as the drift field is 
increased from 1 to 4 kV/cm. This fact has been exploited to obtain a better reconstruction of the 
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deposited energy by combining the two observables. Extensive theoretical work has been devoted 
to understand recombination. Two approaches are typically used. The Jaffe model of recombina- 
tion flZDH assumes that the initial ionization charge is distributed in a column-shaped volume around 
the ionizing particle trajectory. Electrons and ions drift away from this column under the effects 
of drift field and diffusion. During the time evolution of the column, electrons can be captured 
by positive ions, in a process called volume recombination. The Onsager model of recombination 
[UTIl is instead based on the concept of the so-called initial recombination of electron-ion pairs. It 
assumes that the electron has a finite probability of being captured back by the ion Coulomb field, 
and that this probability is reduced by the application of a drift electric field. According to [ITBl , a 
combination of both models is needed to describe the data. 

The transport properties of ionization electrons as they drift in xenon gas, namely the electron 
drift velocity and diffusion, have also been systematically studied as a function of drift field and 
xenon gas pressure (see [E2D and references therein). In the presence of a drift field, electrons 
acquire a net motion in the direction of the field with a stationary drift velocity, v^. Diffusion is 
instead the spread of ionization electrons in the detector volume as the cloud of electrons is drifting. 
These macroscopic properties can be derived by solving the Boltzmann transport equations 112311 . 
provided that a good knowledge of the electron scattering differential cross sections as a function 
of electron energy is available. Different Monte Carlo toolkits for the simulation of electron drift 
and diffusion exist (see, for example, \UM, E5I1 ). To set the scale, for xenon gas at 10 bar pressure 
and for a 1 kV/cm drift field, electron drift velocities of about 1 mm/jUs are obtained, together with 
a longitudinal diffusion spread of a few mm for a 1 m drift. 

No measurable effect on the scintillation signal can typically be attributed to photon transport 
properties in xenon gas. The self-absorption of scintillation photons in pure xenon gas is negligible, 
given that the emission spectrum (due to the formation of electronically excited molecular states, 
or excimers) is well separated from the absorption spectrum of atomic xenon. For this reason, 
any evidence for scintillation light absorption in xenon-based detectors is typically attributed to 
molecular impurities in the gas. 

This paper aims at improving our understanding of some of the xenon gas properties mentioned 
above. Measurements of ionization electron transport properties are described in Sec. |], while 
electron-ion recombination studies are reported in Sec. |5[ This work is performed in the context 
of the NEXT Collaboration R&D, to build a 100 kg scale high-pressure xenon gas time projection 
chamber with the goal of performing a sensitive search neutrinoless double beta decay of the 136 Xe 
isotope flZ5l E7H . Specifically, the data presented in this paper were obtained by exposing NEXT- 
DEMO, a large scale prototype for the NEXT- 100 detector, to a alpha source. Before presenting 
our results (Sec. |] and |5|), we describe the experimental setup in Sec. 0. The NEXT-DEMO data 
processing and selection of alpha candidate events is explained in Sec. |3|. 

2. Experimental setup 
2.1 NEXT-DEMO detector 

NEXT-DEMO is a time projection chamber (TPC) filled with xenon gas at high pressure. The TPC 
fits within a stainless steel pressure vessel, 600 mm long and 300 mm inner diameter, certified for 
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Figure 1. Schematic diagram illustrating the detection principle of the NEXT-DEMO time projection cham- 
ber. 



up to 15 bar operational pressure as well as for operation under vacuum. The detector materials 
place electron-negative impurities into the xenon gas, which need to be removed in order not to 
degrade the detector performance. This is accomplished by continuously re-circulating the xenon 
gas through a getter at a flow of about 10 liters per minute. The TPC fiducial volume consists of a 
hexagonal prism. The hexagon base is defined by PTFE reflector panels, 160 mm across the longest 
diagonal. The fiducial volume extent along the third (drift) dimension, 300 mm long, is defined by 
metallic grids supplying the drift electric field. On one end, the cathode grid is made of 20 stainless 
steel wires at a 10 mm pitch. On the other end, the gate mesh defining the anode is made of 50 
crossed stainless steel wires per inch. The maximum designed drift field is ~ 1 kV/cm. A region 
of electroluminescent (EL, in the following) amplification of the ionization charge is defined by 
creating a high field region between the gate mesh and a second, identical, mesh at ground. The 
passage of ionization electrons in the EL region excite xenon atoms, and a large number of VUV 
photons per ionization electron is produced as a result of excimer de-excitation. The ground mesh 
is placed at a 4.7 mm distance from the gate mesh. The maximum designed EL field between the 
two meshes is ~ 4 kV/(cm- bar). The readout is made of a photomultiplier tube (PMT) plane, 
consisting of 19 Hamamatsu R7378A PMTs, located 120 mm behind the cathode. The PMTs, with 
a 21 mm diameter active area each, are placed in this cathode plane following the hexagonal pattern 
of the reflector panels, 35 mm apart from each other. The active area coverage within the hexagon 
is 40%. These PMTs were chosen because they are capable to withstand 20 bar pressure, and are 
sensitive to VUV light (160-650 nm spectral response). 

As shown in Fig. [I], the PMTs detect both the prompt primary scintillation (SI) light, as well as 
the delayed electroluminescent (S2) light induced by the ionization charge reaching the EL region. 
The time difference between the two light pulses provide the drift coordinate of the event, provided 
that the electron drift velocity is known. The detector can additionally be instrumented with a 
second readout plane, also shown in Fig. [I]. When present, this readout plane is placed a few mm 
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behind the ground mesh, to provide better event position reconstruction in a plane perpendicular to 
the drift direction. Two types of anode readout planes have been tested in NEXT-DEMO: one based 
on 19 Hamamatsu R7378A PMTs, mirroring the cathode readout plane, and another one based on 
~ 250 Hamamatsu S10362-33-050C silicon photomultipliers (see lEBU for details). 

The fast signals produced by the PMTs (less than 5 ns wide) are then shaped and filtered, in 
order to better match the digitizer sampling rate and to to eliminate high-frequency noise. The ana- 
log front-end circuit for NEXT-DEMO is connected via HDMI cables to 12-bit, 40-MHz, digitizer 
cards. These digitizers are read out by FPGA-based data-acquisition (DAQ) modules called front- 
end concentrator (FEC) cards. An additional FEC module with a different plug-in card is used as 
trigger module. Besides forwarding a common clock and commands to all DAQ modules, this card 
receives trigger candidates from the DAQ modules, runs a highly-configurable trigger algorithm 
in the FPGA and distributes a global trigger signal. The trigger electronics accepts also external 
triggers for detector calibration purposes. Whenever a trigger is generated, a PC farm running the 
DAQ software, DATE, receives event data from the DAQ modules via Gigabit Ethernet (GbE) links. 
The DATE PCs (Local Data Concentrators, LDCs) assemble incoming fragments into sub-events, 
which are sent to one or more additional PCs (Global Data Concentrators, GDC). The GDCs build 
complete events and store them to disk for offline analysis. 

2.2 Detector configuration 

All data discussed in this paper were taken during July 23-27, 2012. The detector was filled with 
pure xenon gas at 10 bar pressure. For this NEXT-DEMO run, a flow-through radioactive source 
containing dry radium ( 226 Ra) powder was connected to the gas system. The source, manufactured 
by Pylon (model RNC), has an activity of approximately 74 Bq. The alpha decay of 226 Ra produces 
222 Rn gas that diffuses through the entire detector. As a result, and as shown in Fig. U, three 
main sources of alpha particles are expected within NEXT-DEMO, from the decay of the 222 Rn, 
218 Po and 214 Po isotopes, respectively. The energy deposition of alpha particles is expected to 
be point-like in the detector. More precisely, the CSDA (projected) range of a 5.5 MeV alpha 
particle in xenon at 10 bar pressure (p = 5.56 • 10~ 2 g/cm 3 density) is only 2.4 (2.3) mm. Based 
on previous findings flZQTl . we expect a different spatial distribition of alpha events depending on 
the parent nuclide. Radon is a neutral gas, therefore 222 Rn decays are expected to occur throughout 
the detector. Positively charged 218 Po and 214 Po ions, on the other hand, are expected to (at least 
partially) accumulate on the negatively charged field grids before they are neutralized. 

The NEXT-DEMO detector has been operated with PTFE reflectors both coated with terphenyl- 
butadiene (TPB) and uncoated, see [E8H . The goal of TPB is to shift the xenon VUV light (~ 172 
nm) into blue light (~ 430 nm). This allows for a better light collection efficiency, thanks to the in- 
creased NEXT-DEMO reflectivity as well as the larger PMT spectral response. The data presented 
here were taken with TPB-coated PTFE reflectors. At the time of this run, the PMT-based anode 
readout plane had been dismantled, while the SiPM-based one was still under commissioning and 
not fully operational. For this reason, no anode readout plane information is used in this analysis. 

As for the field configuration, the gate voltage was kept at —5 kV throughout the alpha run, 
resulting in a reduced field between the EL meshes of E^/P ~ 1.06 kV/(cm- bar). This field 
corresponds to a low EL gain of about 200 photons per ionization electron fl3PH . This near-threshold 
EL gain was chosen in order not to exceed the dynamic range of the digitizers, while maintaining 
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Figure 2. Part of the Rn decay scheme that is relevant for this analysis. For alpha decays, if more than 
one energy is allowed, only the dominant one is shown. 

the PMT at their nominal gain of about 5 • 10 6 El- Data were taken at three different drift field 
strengths, Ea^* =0.3, 0.6, and 1 kV/cm, by supplying the cathode with a 14, 23, and 35 kV 
negative high voltage, respectively. 

Data were acquired by triggering on the ionization (S2) signal. The ionization signal charge 
was required to be greater than about 2,000 photoelectrons (PEs) per PMT, in order to reject low- 
energy events. The S2 signal was also required to have an amplitude not exceeding the digitizer 
dynamic range, and a time width of less than 20 /is. The latter requirement ensures keeping all 
alpha events, only a few /is wide, while rejecting tracks with a large (> 2 cm) extent along the drift 
direction, such as most cosmic ray muons. For each trigger, the individual PMT waveforms are 
stored over a 450-600 /is time interval. The DAQ time window was configured in order to have 
100 /is worth of post-trigger data, and pre-trigger data amounting to 350-500 /is, depending on the 
drift field configuration. This time window ensured always recording the S 1 signal corresponding 
to the triggered S2 signal, as it corresponds to more than one full drift length in NEXT-DEMO. The 
average trigger rate was about 5 Hz. Data were acquired for about 24 hours in each of the three 
drift field configurations. 



3. Data processing and selection 
3.1 Waveform treatment 

A pre-processing of the raw data is performed in order to identify the SI -like and S2-like signals 
present in the individual channel waveforms, and to zero-suppress the data. First, pedestals are 
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Figure 3. Typical waveform of an alpha candidate event. This waveform is obtained by summing the 
pedestal-subtracted, gain-corrected, 19 individual PMT waveforms. The SI and S2 signals, appearing also 
in the insets, are clearly visible. 

subtracted for all 19 PMT channels. The pedestal level is computed on a channel-by-channel and 
on a event-by-event basis, using the charge samples within the first 10,000 samples in the waveform 
that are not considered outliers. Second, pedestal- subtracted charge values are converted from ADC 
counts into photo-electrons. These (channel-dependent) gain calibration constants are obtained 
from single-PE LED spectra, from X-ray and photoelectric events from a 22 Na source EBD, as 
well as from alpha events discussed below. Third, the 19 pedestal-subtracted, gain-corrected, PMT 
waveforms are summed together. This summed waveform is less affected by (uncorrelated) noise 
and by statistical fluctuations with respect to the individual channel ones. For this reason, a peak- 
finding algorithm to isolate the SI and S2 signals in the waveform is applied to this channel sum 
only, in a fourth data processing step. 

A typical summed waveform for a alpha candidate event, as defined below in Sec. [Oj, is given 
in Fig. 0. The (early) SI and (late) S2 pulses can be clearly identified. The peak finding algorithm 
first builds pulses, or "peaks", as collections of consecutive samples above pedestal. Peaks can 
qualify as Sl-like or S2-like depending on their time and charge characteristics. A Sl-like peak is 
required to be less than 3 /J.S wide, and to have an average charge per PMT of at least 0.5 PEs/PMT. 
A S2-like peak is defined as being at least 3 fis wide, and with a charge of at least 500 PEs/PMT. 
Because of their wide pulse structure, for S2 peaks only the charge samples are added together 
into 0.25 jis time bins (that is, by applying a resampling factor of 10). This allows for some data 
reduction at no significant information loss. Fifth, the bulk of the raw data reduction is obtained by 
dropping from the individual channel waveforms all samples that are not coincident in time with 
the Sl-like and S2-like peaks. 

3.2 Selection of alpha-like ionization signals 

The first step in the event selection, common to both alpha samples discussed in the following 
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Table 1. Data reduction summary for the selection of alpha-like ionization signals. The number of events 
before and after each selection requirement is given, for all three drift field configurations studied. 



Requirement 


Drift Field (kV/cm) 
0.3 0.6 1.0 


None 
1 S2 
1 alpha-like S2 


361,796 488,017 397,391 
347,103 478,203 368,974 
341,981 431,052 329,226 



(Sees. ^3] and [3~31 ), involves the analysis of the S2 signal only. We first require a single S2-like 
peak (as defined in Sec. [3~T1 ) in the entire channel sum waveform of the event. This requirement 
is motivated by suppressing complicated topologies, such as the ones involving accidental coinci- 
dences. We then require the S2 peak to be alpha-like in order to select the event. A S2 peak is 
defined as alpha-like if it fulfills two conditions. First, the charge detected in time coincidence with 
the S2 peak must exceed 1000 PEs for all of the 19 PMTs in the plane. This is to ensure that the 
S2 peak on the channel sum waveform is a genuine S2 signal, illuminating all PMTs at once, as 
opposed to uncorrelated noise on only some of the PMTs. Second, reinforcing what already speci- 
fied at the trigger level, the S2 peak is required to have a narrow time width comprised between 2 
and 18 jJ.s, where the width is defined in this case as the full width at 5% of the S2 pulse height. 
In other words, only energy depositions localized in 3D space (such as alpha decays) or extending 
only perpendicularly to the drift direction are expected to fulfill this requirement. A data reduction 
summary for the selection of alpha-like ionization signals is given in Tab. [I]. 

3.3 Selection of alpha decays in the gas 

In addition to have alpha-like ionization signals, alpha candidate events in the xenon gas are re- 
quired to fulfill the following criteria on the scintillation (SI) signal. First, at least one SI -like 
peak (as defined in sec. [3~T1 ) is required on the channel sum waveform. At this stage we allow for 
more than one SI peak to avoid large selection efficiency losses, given that multiple spurious SI 
peaks are often found by our peak finding algorithm. Second, we require exactly one pulse, among 
this pool of S 1 candidate pulses per event, to be alpha-like. A alpha-like S 1 signal is defined as 
one having a signal risetime of less than 200 ns, and a minimum charge of 20 PEs per PMT in 
time coincidence with the SI peak, for all 19 PMTs in the readout plane. The risetime requirement 
eliminates S2-like pulse shapes, which are more symmetric and with longer risetimes compared to 
SI -like pulse shapes (see insets in Fig. |3|). The PMT multiplicity requirement ensures that genuine 
SI signals are selected as opposed to noise, as well as sufficiently bright ones. 

At this stage, the combination of the event SI and S2 signals allows us to reconstruct the alpha 
decay candidate in 3D space within the chamber. Our coordinate system is defined as follows: z 
axis along the drift direction and pointing from anode toward cathode, y axis pointing upward, x 
axis forming a right-handed coordinate system. The origin is located at the intersection between 
the chamber axis and the plane parallel to the EL meshes and located in between them. 

As a third condition, we require the event drift length z to be comprised between 40 and 220 
mm, as shown in Fig. |J. The drift length is computed by multiplying the event drift time by the 
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Figure 4. Drift length distribution of events in the i^jrift = 0.6 kV/cm run, prior to the drift length cut. 
Events between the two vertical lines at 40 and 220 mm drift are retained. The large spike at ~ 300 mm drift 
corresponds to events originating from the cathode plane. 



drift velocity measured for each drift field configuration, as discussed in Sec. fOl . In turn, the drift 
time is given by the difference between the SI ans S2 peak positions. As can be seen in Fig. |], 
alpha decays occur throughout the full, 300 mm long, drift volume l . The large accumulation of 
events near 300 mm corresponds to alpha decays from the cathode plane, which we discuss in 
Sec. PI Events at drift distances smaller than 40 mm are difficult to reconstruct because of the 
proximity of the SI and S2 signals. Events at drift distances larger than 220 mm, and up to 300 
mm, can instead be reconstructed. However, they are excluded from this analysis because they 
involve a more complicated S 1 yield spatial correction with respect to the linear one applied here, 
see Sec. Q 

The fourth and final condition retains events whose reconstructed (x,y) position is within a 
central, 30 mm radius, cylindrical region. For this run, reconstruction in the (x,y) plane is com- 
plicated by the lack of a readout plane near the EL region. The isotropically emitted EL light will 
in fact illuminate in a almost uniform manner the cathode plane readout, located 420 mm away. 
However, minimal (x,y) position information is retained on this plane, for very bright events such 
as alpha decays, and provided that the PMT responses are properly inter-calibrated. More quanti- 
tatively, from solid angle arguments only and assuming perfectly absorptive reflector panels, PMT 
response inhomogeneities of the order of 35 2 /420 2 ~ 0.7% are to be expected between a PMT at 
the same (x,y) location as the light production point in the EL region, and a neighbouring PMT 
located 35 mm away. This number is to be compared with the expected level of statistical fluctua- 
tions. As we will see in Sec. [J/j] , alpha candidate events produce S2 signals with a typical charge 
of 5,500 PEs per PMT. Therefore, while sub-dominant, inhomogeneities due to solid angle angle 
effects are comparable with statistical fluctuations, at the (5500) -1 / 2 ~ 1.3% level. We define the 

1 The depletion of events at very small drift distances, z < 20 mm, is due to an artifact of the peak finding algorithm. 
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Figure 5. Reconstructed position of events in the (x,y) plane for the i^jrift = 0.6 kV/cm run, prior to the 
radial cut. Events within the central, 30 mm radius, region are retained. An accumulation of events at the 
position of the 19 PMTs can be seen, and is due to a reconstruction bias. 

event (x,y) position via a charge-weighted sum of the PMT positions, using a modified barycenter 
(or centroid) method: 

N 

* = E((<7//4)7A0-^ (3-1) 

i=\ 

where N = 19 is the total number of PMTs, q is the average S2 charge for that event, Xi is the x 
position of PMT i, and qt is the S2 charge measured by PMT i (and similarly for the y coordinate). 
In the normal barycenter algorithm, the exponent n is set equal to 1 . However, the near-uniform 
illumination of the PMT plane would largely bias the (x,y) reconstruction toward (0,0), in this 
case. For this reason, we modify the method to greatly amplify inhomogeneities at the level ex- 
pected by solid angle effects, by choosing a large value for the exponent in Eq. PTT] , n = 200. This 
value was guided by the requirement 1.007" ~ 4, where 4 is approximately the average number 
of neighbouring PMTs in the readout plane. Events satisfying the requirement r = ^Jx 1 +y 2 < 30 
mm, where x and y are defined as in Eq. PTTj , are retained in our final sample of alpha candidate 
events in the gas, as shown in Fig. [5]. 

A summary of the data reduction, from alpha-like ionization signal events (Sec. [3T2D to alpha 
candidate events in the gas, is given in Tab. 0. 
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Table 2. Data reduction summary for the selection of alpha decays in the gas. The number of events 
surviving each selection requirement is given, for all three drift field configurations studied. 



Requirement 


Drift Field (kV/cm) 




0.3 


0.6 


1.0 


1 alpha-like S2 


341,981 


431,052 


329,226 


>1 SI 


287,797 


357,797 


262,916 


1 alpha-like SI 


263,809 


334,041 


241,444 


Drift length cut 


156,950 


184,043 


134,926 


Radial cut 


58,737 


75,162 


57,834 




Drift Length (mm) 

Figure 6. Dependence of the S 1 signal yield as a function of drift length. The 2D histogram of gas alpha 
candidate events for the i^jrift = ^-6 kV/cm run (colored boxes), the corresponding profile histogram (black 
points), and a linear fit to the profile histogram (curve) are shown. 

3.4 Spatial corrections to ionization and scintillation signals 

The ionization and, especially, the scintillation yields are affected by well-understood instrumental 
effects, causing the observed yields to depend on the alpha decay spatial position within the gas. 
The most important, and easiest to correct, effect is given by the dependence of both the ionization 
and scintillation signals with drift length. This correction is discussed below, using the sample 
of gas alpha candidate events (see Sec. P3| ). There is also a radial dependence of the yields, as 
discussed in EBD- Given the difficulty in reconstructing the radial position of the decays without 
a anode plane readout, for the data discussed here we do not attempt any radial correction. As 
mentioned above, only a simple radial cut using cathode plane information only is applied to our 
data, in order to retain mostly central events. 

The dependence of the scintillation signal with drift length is shown in Fig. As the alpha 
decay occurs closer to the energy plane located behind the cathode (large drift lengths), more 
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Figure 7. Dependence of the S2 signal yield as a function of drift length. The 2D histogram of gas alpha 
candidate events for the i^jrift = 0.6 kV/cm run (colored boxes), the corresponding profile histogram (black 
points), and a exponential fit to the profile histogram (curve) are shown. 

scintillation light is detected on average. The reason is the non-perfect reflectivity of the reflector 
panels. In other words, while the TPB coating greatly improves the light collection efficiency and 
reduces its spatial dependence, it does not eliminate this dependence completely. In the fiducial 
volume region chosen for this analysis, between 40 and 220 mm drift length, the average S 1 yield 
increases by about a factor of 2, increasing approximately linearly with drift length. The S 1 yields 
shown in the following (Sec. 0) are corrected according to this linear dependence, to provide a 
S 1 yield at zero drift length. While the data in Fig. ^| refer to 0.6 kV/cm drift field data only, we 
have verified that this spatial dependence is independent of drift field strength, as expected from 
the above considerations. 

As shown in Fig. [7], also the average ionization yield varies with drift length. In this case, the 
dependence is opposite, with more detected charge closer to the EL region (short drift lengths). 
Impurities in the xenon gas, causing some amount of attachment of ionization electrons, are re- 
sponsible for this effect. The 0.6 kV/cm data shown in Fig. |7J are well described by an exponential 
function, corresponding to a 1.2 ms electron lifetime, and to a ~15% reduction in the average S2 
yield between 40 and 220 mm drift length. Unlike for the S 1 yield, this correction is dependent on 
the drift field configuration. This attachment effect is the largest in the lowest drift field configura- 
tion (0.3 kV/cm), where the average S2 variation across the 40-220 mm drift length range becomes 
~20%. The S2 yields shown in Sec. [5] are corrected according to this exponential dependence, to 
provide a S2 yield at zero drift length. 

3.5 Selection of alpha decays on cathode wires 

As shown in Fig. |], several alpha candidate events originate from the cathode plane at z — 300 mm. 
A clean sample of cathode alpha candidate events can be easily isolated from alphas decays the gas, 
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Table 3. Data reduction summary for the selection of alpha decays in the cathode plane. The number of 
events surviving each selection requirement is given, for all three drift field configurations studied. 



Requirement 


Drift Field (kV/cm) 
0.3 0.6 1.0 


1 alpha-like S2 

>1 SI 
1 cathode-like S 1 


341,981 431,052 329,226 
287,797 357,797 262,916 
25,829 50,363 33,327 



starting from the alpha-like ionization signal selection (Sec. [3T2D . As in Sec. [33j we first require 
at least one SI -like peak. Among those peaks, we retain only the ones with a < 200 ns risetime 
(as in Sec. [3"3D and an associated charge on each of the 19 PMTs of at least 200 PEs. We define 
as cathode alpha candidates the events with one and only one S 1 peak fulfilling the two conditions 
above. These events are therefore characterized by abnormally large SI signals, above the level 
expected for SI signals from alpha decays in the gas, see Fig. ^. In addition, we have observed that 
the average SI yield for the cathode alpha candidate events increase very markedly with the high 
voltage supplied to the cathode, from ~ 300 PEs/PMT at -14 kV to ~ 1, 500 PEs/PMT at -35 kV. 
These observations can be explained by having some amount of electroluminescence produced on 
the cathode wires. On the other hand, the S2 yield for such events is similar to the one for alpha 
particles in the gas. In this case, the S2 charge distribution exhibits a bimodal distribution which 
is consistent with being due to 218 Po (6.00 MeV) and 214 Po (7.69 MeV) alpha decays. An event 
selection summary for cathode alpha candidate events is given in Tab. || 

4. Measurement of ionization electron transport properties 
4.1 Electron drift velocity 

The sample of alpha decays from the cathode wires (cathode alpha candidate events, in the follow- 
ing) allows for a very accurate determination of the electron drift velocity. A drift time distribution 
for such events is shown in Fig. |8|. Specifically, the events shown in the figure refer to the first 
two hours of data taken in the = 0.6 kV/cm configuration. The data are well described by 

a gaussian distribution (also shown in Fig. ||) with a mean drift time of (299.848 ±0.023) jJ.s and 
a one sigma spread of 0.9 [is. The spread in drift times for these events is therefore very small. 
Possible factors contributing to this spread include noise, the relatively wide (few jits) ionization 
signal time extent, the orientation of the alpha decays with respect to the drift direction, sub-mm 
non-uniformities in the NEXT-DEMO full drift distance as a function of radial position, and drift 
velocity variations with time (see below). 

We obtain a measurement of the drift velocity from the ratio between a fixed drift distance 
for all cathode alpha candidate events, and a mean drift time as the one shown in Fig. |8[ The 
relevant distances are the 300 mm distance between the cathode and the gate grids, the 4.7 mm 
distance between the gate and the ground grids {eg, the width of the EL region), and the 1.15 
mm average range projected along the drift direction for cathode alpha candidate events emitted 
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Figure 8. Drift time for cathode alpha candidate events. Only data from approximately the first two hours 
of the i^jrift = 0.6 kV/cm run are shown. The red curve is a gaussian fit to the data. 
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Figure 9. Drift velocity as a function of time, for the three drift field configurations studied. For each time 
interval, the data points and errors shown are obtained from a gaussian fit to the drift time data as the one 
in Fig. H, plus the known full drift length of the NEXT-DEMO detector. The three small panels on the right 
show a zoomed version of the same drift velocity data as the main panel. 



in the hemisphere enclosing the detector fiducial volume. Recalling that the drift time is given 
by the ionization signal peak time minus the scintillation signal peak time, and assuming that the 
ionization signal peak occurs in the middle of the EL region, we take 300 + 2.35 — 1.15 = 301.2 
mm for this fixed drift distance. 

The variation in the drift velocity measurements as a function of time within a given run, and 
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Figure 10. Drift velocity as a function of drift field. The black points are the measured values, the solid 
(dashed) curve is the prediction for pure xenon from the Magboltz \ UM (Escada et al. l~ E5ll ) simulation. The 
Magboltz results are based on version 9.0. 1 of the code, and include thermal motion of atoms. The error bars 
on the NEXT-DEMO data account for the main uncertainties (drift length and pressure, see text for details). 

for approximately two hour time intervals, is shown in Fig. |9| for the three drift field configurations 
studied. The main aspect to highlight is the excellent stability over time, with only <0.5% drift 
velocity variations during more than one day of data. Even so, the accuracy of the measurements 
is sufficient to appreciate a systematic trend in the data, in the form of a very slow drift velocity 
increase over time, at least for the 0.6 and 1.0 kV/cm runs. This small increase in drift velocity 
over time (on average) can be explained by a gas leak that was observed (by manual inspection 
of the pressure gauges) during the run. The observed pressure leakage average rate was about 
0.3%/day. Gas was routinely (every few days) added to the chamber to compensate for the gas 
leak but, unfortunately, no detailed track record of the pressure conditions during the run exists. 
From the above considerations, in the following we will assume a 2% systematic uncertainty on 
the pressure conditions in NEXT-DEMO during this run. 

We have also compared our measured drift velocities at 0.3, 0.6, and 1 kV/cm drift fields with 
expectations from pure xenon gas at 10 bar pressure obtained with two independent simulations of 
electron transport in gases [E3L E3- The comparison is shown in Fig. |TtJ- The agreement is fair 
for all drift field configurations and for both simulations, with observed drift velocities only a few 
% larger than the expected ones for pure xenon. The observed drift velocity central values were 
taken by simply averaging over the entire run at a given drift field. The error bars shown in Fig. [K^ 
assume two main contributions, summed in quadrature: a 2% uncertainty in the NEXT-DEMO full 
drift length (directly translating into a 2% drift velocity uncertainty), and the 2% pressure reading 
uncertainty mentioned above (translating into a ~1% drift velocity uncertainty). 

4.2 Electron longitudinal diffusion 

The diffusion of ionization electrons as a function of drift distance can instead be studied with alpha 
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Figure 11. Dependence of the S2 pulse time width (FWHM r ) as a function of drift length, illustrating the 
effect of longitudinal diffusion. The 2D histogram of gas alpha candidate events for the ^rift = ^-3 kV/cm 
run (colored boxes), the corresponding profile histogram (black points), and a fit to the profile histogram 
according to Eq. [CT] (curve) are shown. 



decays in the xenon gas (gas alpha candidate events, in the following), by exploiting the point-like 
nature of alpha particle energy depositions within the chamber. As mentioned above, the range of 
alpha particles is expected to be of the order of 2 mm, roughly of the same order of (and often 
smaller than) diffusion effects. 

In this work, we infer the longitudinal diffusion of ionization electrons in xenon gas by study- 
ing how the ionization signal pulse width stretches with increasing drift length. In our simple 
model, we fit the full width, half maximum of the S2 pulse time extent (FWHM r ) via the following 
function: 



FWHM, = 2V21n2 • J o% + of (4. 1) 

where do is a constant term, and the longitudinal diffusion term Of,, which varies with drift length, 
is defined as: 



al = (2D L /v 3 d )-z (4.2) 

where z is the drift length, v c \ is the electron drift velocity, and Di is the longitudinal diffusion 
coefficient (in cm 2 /s units). In our model, there are two free fit parameters, o and Di, while vj is 
taken from our measured values using cathode alpha candidate events (Sec. jOp . 

Figure |TT] shows the dependence of FWHM, in Eq. ^71] with drift length, z. For illustrative 
purposes, we show here the data for the 0.3 kV/cm drift field configuration, that is our lowest field 
measurement, in order to maximize the effect of diffusion. The 2D histogram of FWHM, versus 
drift length for all gas alpha candidate events shows a clearly defined band, indicating a very pure 
sample of correctly reconstructed alpha candidate events. In other words, few mis-reconstructed 
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Figure 12. Longitudinal diffusion coefficient Dl as a function of drift field. The black points are the 
measured values, the solid (dashed) curve is the prediction for pure xenon from the Magboltz tUM (Escada 
et al. r E5ll ) simulation. The Magboltz results are based on version 9.0.1 of the code, and include thermal 
motion of atoms. The error bars on the NEXT-DEMO data account for the main uncertainties (drift velocity 
and fit model, see text for details). 

events due to accidental S1+S2 combinations, and few events with extended energy depositions 
along the drift direction, survive our event selection described in Sec. P3| . Moreover, the width of 
this band, of order 1 pLs for all drift lengths, is in qualitative agreement with expectations for the 
non-zero alpha range in xenon gas at 10 bar pressure. 

The profile histogram of FWHM, versus drift length, also shown in Fig. [TT], is used as input to 
the fit given by Eq. fPl . Our simple model, whose best-fit line is shown in Fig. [TT] as well, provides 
a fair but far from perfect description of the data, with % 2 values per degree of freedom of about 
10. 

Figure [J2] shows how our best-fit longitudinal diffusion coefficients Di for the 0.3, 0.6, and 
1 kV/cm field configurations compare to the predicted values for pure xenon gas at 10 bar, ob- 
tained from the same two simulations discussed in Sec. |Q| . The spread in the two predictions is 
much larger than the corresponding one for the drift velocity predictions (Fig. UJJ). Our data are in 
qualitative agreement with the prediction for pure xenon from Escada et al. IE511 . clearly favoring 
the latter with respect to the Magboltz (version 9.0.1, atomic thermal motion included) prediction 
[EH]. The error bars on the NEXT-DEMO diffusion data assume two main contributions, summed 
in quadrature: a drift velocity uncertainty as estimated in Fig. [TTJ, plus a fit model uncertainty, 
obtained by scaling the fit uncertainties by y^/dof. 

Finally, the best-fit constant term a in Eq. ^TJJ shows little dependence with drift field, as 
expected, and values comprised between 0.80 and 0.93 pts for the three field configurations studied 
are returned by the fit. Factors contributing to a non-zero ionization signal pulse time extent at 
zero drift length (or, equivalently, for negligible diffusion) include the shaping in time of the signal 
due to the electronics, the 4.7 mm spatial extent of the EL region, and the non-zero alpha range in 



xenon gas at 10 bar pressure. 

5. Ionization yield, scintillation yield, and electron-ion recombination studies 

In Sec. |], we use observables in the time domain (namely, drift time and time extent of the ioniza- 
tion signal) to infer ionization electron transport properties. In this section, we use instead charge 
information to study the scintillation yields, the ionization yields, and the relationship between the 
two. As for the diffusion measurements in Sec. fOj we use the sample of gas alpha candidate 
events described in Sec. E3- 

5.1 Correlated fluctuations between ionization and scintillation 

Figure [T3| shows the ionization (S2) yield as a function of the scintillation (SI) yield, for gas alpha 
candidate events collected in the 0.6 kV/cm field configuration. Both yields have been corrected 
for the (detector-dependent and impurity-dependent) spatial effects described in Sec. ^3]. The pro- 
jections along the two axes are also shown in Fig. [T3[ As hinted by the 2D histogram (and more 
clearly demonstrated below), the majority of the events correspond to alpha decays of 222 Rn (5.49 
MeV energy, see Fig. §), plus a ~ 15% fraction of alpha decays of 218 Po (6.00 MeV energy). Large 
fluctuations are visible both in the ionization and scintillation response, spoiling the identification 
of the 222 Rn and 218 Po alpha lines using the SI or S2 signals alone. Despite the fact that the number 
of photo-electrons detected for the ionization signal is about 50 times larger than the corresponding 
number for scintillation (see Figs. |5| and (7)), fluctuations are larger for the ionization channel. How- 
ever, the data also show a clear anti-correlation between the ionization and scintillation signals. 
For this reason, we adopt a similar procedure as the one followed in liquid xenon experiments to 
combine the two signals, in order to obtain a more accurate estimator of the deposited energy per 
event. Our procedure, described below, closely follows the one in HT9H . 

In order to estimate the alpha energy, we proceed in the following, iterative, way. In a first step, 
we use the mean SI and S2 values of the 2D histogram in Fig. [T3], in PEs/PMT units for the SI 
and S2 corrected yield at zero drift length, to rescale the SI and S2 response into absolute (MeV) 
energy units. In the following, we refer to these MeV-rescaled SI and S2 yields as E$\ and E$i, 
respectively. At this stage, we assume that all events are 5.49 MeV alphas from 222 Rn decays. In a 
second step, we fit the entire 2D (£"51 , £52) histogram shown in Fig. [13] to a 2D correlated gaussian 
distribution. The fit therefore includes six parameters: an overall normalization, mean and sigma 
of the SI signal in MeV units, mean and sigma of the S2 signal in MeV units, and a dimensionless 
correlation coefficient p between the two. In a third step, we estimate the S1-S2 correlation angle 
# as: 

■& = arctan (052/051 ) (5.1) 

where 051 and 052 are the gaussian fit sigma parameters along the two axes. In other words, ■& is 
the angle with respect to the S 1 axis of the major axis of the ellipse obtained from the 2D gaussian 
fit. The red solid line in Fig. |T3| is defined as the line having this slope ■& with respect to the SI 
axis, and passing through the mean £"51 and £52 parameters returned by the fit. In a fourth step, we 
estimate the event energy E by performing a rotation in the (E$\ ,£52) plane according to this angle 
#, as follows: 
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Figure 13. Ionization (S2) signal yield as a function of scintillation (SI) signal yield for gas- alpha candidate 
events in the ^rift = ^-6 kV/cm run (colored boxes), illustrating the effect of correlated fluctuations between 
them. The two yields have been corrected for the spatial effects described in Sec. [3.41 , and also normalized 
to an absolute energy (MeV) scale. The events around Eg^ = 5.5 and 6.0 MeV are mostly Rn and Po 
alpha decay events, respectively. The red curve is the E = 5.49 MeV line, where the energy E is defined as in 
Eq. |5.2| . In order to extract the red line, only events to the left and bottom of the dashed (E = 6.0 MeV) line 
are used. Finally, the ID projections on the ionization and scintillation axes are also shown as histograms. 



E = : — (5.2) 

sint> +cos n 

In a fifth step, we apply an energy cut to partially remove the biases introduced by the fact 
that not all events in our sample are alpha decays from 222 Rn. In order to do so, we repeat the 
steps above only for events satisfying the requirement E < 6 MeV, where E is defined as in Eq. |5T2"1 . 
This is the black dashed line shown in Fig. |T3|. The 6 MeV value for this cut was chosen as a 
compromise between maintaining a high efficiency for 222 Rn alphas while rejecting as many 218 Po 
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Table 4. Values of the parameters Gsi , G52 and p for the 2D gaussian fit to gas alpha candidate events in 
the (Esi,Es2) plane, for the three drift field configurations studied. The value of the angle # as defined in 
Eq. |5TT] is also given. 



Drift field 




&S2 


P 


1? 


(kV/cm) 


(MeV) 


(MeV) 




(rad) 


0.3 


0.307 ±0.001 


0.595 ±0.003 


-0.291 ±0.006 


1.094 ±0.001 


0.6 


0.251 ±0.001 


0.574 ±0.003 


-0.221 ±0.006 


1.158 ±0.001 


1.0 


0.234 ±0.001 


0.579 ±0.003 


-0.201 ±0.007 


1.186±0.001 




Figure 14. Energy distribution for gas alpha candidate events in the = 0.6 kV/cm run, where the 

energy is defined according to Eq. |5.2| . Two peaks, due to 222 Rn and 218 Po alpha decay events, are visible. 
The red curve is a fit to the data, assuming the sum of two gaussians to describe the two populations. 



events as possible (about half). One or two additional iterations of this entire procedure have proven 
to be sufficient to obtain stable numerical results. 

The parameters Ogi* Gs2 and p of the 2D gaussian fit, together with the derived quantity ■&, are 
reported in Tab. |] for the three drift field configurations studied. For all drift fields, the correlation 
coefficient is negative, and the resolution using the SI signal alone is about a factor of 2 better 
than using the S2 signal alone. As expected, and as already observed in liquid xenon [TT5L 033], the 
anti-correlation is strongest for the lowest drift field configuration, corresponding to the strongest 
recombination scenario. 

The event energy E for all gas alpha candidate events collected during the 0.6 kV/cm drift 
field run is shown in Fig. [l4]. Two peaks are clearly visible, corresponding to 222 Rn and 218 Po 
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Table 5. Values of the parameters for the double gaussian fit to the energy spectra, for the three drift field 
configurations studied. The mean and sigma of the 222 Rn ( 218 Po) gaussian are indicated with Ep^ and Op n 
(£p and 0p o ), respectively. The ratio of 218 Po-to- 222 Rn events is given by Np /Np n . 



Drift field 


%n 


a Rn 


^Po 


a Po 




(kV/cm) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 




0.3 


5.476 ±0.002 


0.198 ±0.004 


6.034 ±0.010 


0.193 ±0.005 


0.208 ±0.010 


0.6 


5.483 ±0.002 


0.1 80 ±0.002 


6.036 ±0.006 


0.1 87 ±0.004 


0.211 ±0.006 


1.0 


5.488 ±0.002 


0.170 ±0.002 


6.025 ±0.007 


0.193 ±0.004 


0.209 ±0.006 



alpha decays, respectively. For energies above 5.3 MeV, the data are well described by a sum of 
two gaussians for these two populations. This fit is also shown in Fig. 14. While the 222 Rn peak 



position is essentially required by our procedure to be aligned with its expected value of 5.49 MeV, 
this is not the case for the 218 Po peak. The fact that the latter appears near the expected 6.00 
MeV value gives us confidence in our hypothesis that most events in our sample are 222 Rn alpha 
decays, with about a 15% fraction due to 218 Po alpha decays. Both gaussian distributions in the fit 
of Fig. [T4| are characterized by a FWHM energy resolution of about 8%, slightly improving with 
increasing drift field, about a factor of 3-4 worse than what has already been obtained with the 
same NEXT-DEMO detector (albeit in a somewhat different configuration) using a 22 Na 511 keV 
gamma ray source, see [12811 . 

As can be seen from the fit parameters listed in Tab. |5], the energy spectra do not change in 
any appreciable manner as the drift field is increased from 0.3 to 1.0 kV/cm. As mentioned above, 
the depletion of 218 Po decays in the gas with respect to 222 Rn decays seems to indicate that the 
majority of the positively charged 218 Po ions accumulate on the negatively charged field grids. The 
214 Po decays are even further suppressed in our gas alpha sample selection, since the very short 
half-life of their parent 214 Bi isotopes (see Fig. ||) implies that our requirement of a single SI and a 
single S2 signal per event (see Sees. and [331 ) is rarely fulfilled in this case. 

The energy distribution in Fig. [14] also shows a non-gaussian low-energy tail, extending down 
to 4.5 MeV energies and beyond. This low-energy tail is caused by the events with abnormally low 
(E S 2 ~2-4 MeV) S2 yields in Fig. [T3[ The origin of this tail, present for all drift field configurations, 
is not understood at the moment. Early investigations seem to exclude ionization charge loss on the 
NEXT-DEMO reflector walls as a cause. This tail is in fact only mildly affected by requiring and/or 
tightening the radial cut selecting central events, both for alpha decays in the gas discussed here as 
well as for alpha decays in the cathode plane 2 . However, a more definitive test of this hypothesis 
will be made once additional data will be taken, operating NEXT-DEMO with an additional anode 
readout plane. 



2 The investigation of cathode alpha events is interesting in this respect, given that a better (x,y) reconstruction 
compared to gas alpha events can in principle be obtained in this case. This was accomplished by using the large SI 
signals (as opposed to S2 signals) produced at a fixed distance from the PMTs, emitted 120 mm (as opposed to 420 mm) 
away from the PMTs. 
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Figure 15. Stability over time for the average S 1 and S2 signals of Rn events, for the three drift field 
configurations studied. The filled black data points assume a spatial correction to the signals (see Sec. |3.4| ) 
that is constant throughout a drift field configuration. The empty data points for the S2 stability over time 
assume a time-dependent spatial correction. 

5.2 Field dependence of ionization and scintillation average yields 

Apart from the study of correlated fluctuations discussed in Sec. [5TTT, electron-ion recombination 
in xenon gas can also be investigated by measuring how the average ionization and scintillation 
yields vary with drift electric field intensity. For very intense fields, recombination is negligible. 
In this case, ionization charge collection is unaffected by the charge carriers in the gas, and the 
scintillation signal is entirely due to primary scintillation. As the field is lowered and recombination 
becomes important, part of the ionization is lost while the scintillation signal increases because of 
recombination luminescence. Different recombination models yield different trends in the way 
both signals vary with electric field intensity and with ionization density (see, for example, HTBl ). 
For the drift field and gas pressure conditions probed here, variations at the few % level in both the 
average SI and S2 yields are expected, see for example UTTH . 

The necessity for such accurate (% level) studies of S 1 and S2 variations with drift field require 
instrumental effects to be understood at the same level or better. Figure |T3| shows the stability over 
time of the average S 1 and S2 yield for 222 Rn alpha events, for all three drift field configurations 
probed. The yields has been corrected for spatial effects, as explained in Sec. [ITtl More specifically, 
the average yields are obtained as the center of the ellipse in the {Esi,Es%) plane obtained from 
the 2D gaussian fit to the 222 Rn alpha decay events, as explained in Sec. However, unlike 
Fig. |T3|, the yields have not been rescaled into absolute energy units. The S 1 response is stable to 
within 0.4% or better, for all drift field configurations. On the other hand, the S2 response tends 
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Figure 16. Average scintillation yield as a function of drift field, for gas alpha candidate events. 



to increase with time. This is consistent with pressure slightly decreasing over time because of a 
gas leak in the system, as discussed in Sec. |7T]. In fact, the S2 yield trends over time in Fig. |T3| are 
similar to the drift velocity trends in Fig. 0, with a more (less) evident slope for the 1.0 kV/cm (0.3 
kV/cm) data. Compared to drift velocity, however, the S2 yield is expected to be more sensitive 
to small pressure variations. More quantitatively, a 2% pressure decrease at 10 bar (the size of 
our estimated uncertainty on the pressure reading) is expected to generate about a 4% increase in 
the S2 response [BQl . We can in principle correct for such pressure variations during the data- 
taking period by assuming a time-dependent spatial correction to the S2 signal, albeit with larger 
stattical fluctuations compared to a time-independent correction, as shown in Fig. [T3|. However, 
and more importantly, pressure variations of the same order may have happened while the detector 
was not taking data. For this reason, while we consider the S 1 response to be understood at the 
level necessary for this study, the S2 response is not. 

Figure 16 shows how the average scintillation yield for the 222 Rn line varies with drift field. 
The points and error bars in the figure represent the mean and the standard deviation of all measured 
time intervals at a given drift field in Fig. [13]. We observe that the average scintillation yield 
increases by (10.6±0.5)% as the field is decreased from 1.0 to 0.3 kV/cm. This variation appears 
to be similar (perhaps somewhat larger) to what found by the authors of [11711 for xenon gas in 
the same pressure conditions, and also using alpha particles. A similar study will be repeated in 
the future also for the ionization signal, by operating NEXT-DEMO in more stable (and better 
monitored) conditions. 



6. Conclusions 

During the past few years, the NEXT Collaboration has conducted an R&D program geared to- 
ward the construction of a high-pressure xenon gas time projection chamber (TPC) at the 100 kg 
scale (NEXT- 100), to perform a sensitive search for neutrinoless double beta decay in 136 Xe. Sev- 
eral time projection chamber prototypes have been built in this context, to prove the technological 
viability of the proposed NEXT- 100 detector concept and to reliably quantify its projected perfor- 
mance. NEXT-DEMO is the largest of such prototypes, and has proven to be an ideal tool to also 
study the detection properties of high-pressure xenon gas in a broader context. 
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In this paper, we have studied both the ionization and the scintillation response of 10 bar pres- 
sure xenon gas to alpha particles, operating the NEXT-DEMO time projection chamber at 0.3-1.0 
kV/cm drift field strengths. Alpha particles from the decay of the 222 Rn, 218 Po and 214 Po isotopes 
are introduced in the detector, by allowing the xenon gas to flow through a 226 Ra radioactive source. 
Two clean samples of alpha candidate events are used, one from decays in the xenon gas and one 
from decays from the TPC cathode, operated at a large negative high voltage. 

Alpha particles are used to study both the transport properties of ionization electrons in the gas 
as well as the process of electron-ion recombination in xenon gas. Concerning electron transport, 
an accurate monitoring of electron drift velocity was obtained using alpha decays from the cathode 
surface, confirming model predictions for pure xenon to the few percent level. The longitudinal 
diffusion of ionization electrons as they drift has also been accurately measured, by exploiting the 
point-like nature of alpha-induced energy depositions in the gas. Longitudinal diffusion is less well 
understood than electron drift velocity, and different transport models yield longitudinal diffusion 
coefficients, differing by as much as a factor of two imil . Our data clearly favor low diffusion 
models, such as [C5I1 . This result has a direct impact on the NEXT- 100 physics potential, since the 
signal-to-background discrimination achievable in the detector by using event topology is affected 
at some level by electron diffusion. 

The process of electron-ion recombination can be best understood by the simultaneous mea- 
surement of both the scintillation and the ionization response, as in NEXT-DEMO. When recom- 
bination is negligible, the observed scintillation light is entirely primary, and no ionization charge 
losses (apart from attachment) occur. As electron-ion recombination becomes non-negligible, part 
of the ionization signal is lost to the benefit of a stronger scintillation response. The amount of 
recombination is expected to increase with decreasing electric field and with increasing linear en- 
ergy transfer, the latter being defined as the quotient between the particle energy deposition per unit 
path length and the density of the detection medium. Despite the relatively low density of 10 bar 
pressure xenon gas (that is, low density relative to liquid xenon), the use of highly ionizing alpha 
particles allows us to clearly observe recombination in xenon gas in two different ways. On the 
one hand, we observe a quenching of the average scintillation response of xenon gas as the field is 
increased, as already observed several times in the past f fl"5l . ITfil H7I] . The scintillation light response 
to 5.49 MeV alpha particles from 222 Rn decreases by about 10% by increasing the drift field from 
0.3 to 1.0 kV/cm. A similar result has been obtained in [H7H . On the other hand, and perhaps more 
interestingly, we observe large event-by-event correlated fluctuations between the ionization and 
scintillation response. We measure negative correlation coefficients, and a stronger anti-correlation 
for lower drift fields, as expected from recombination. We exploit this anti-correlation to build a 
better energy estimator that combines both the ionization and scintillation response. While similar 
findings and strategies have been reported in the past in liquid xenon fTTEl . 033], to our knowledge 
this is the first time that event-by-event correlated fluctuations between ionization and scintillation 
are observed in xenon gas. 

In addition to better understand the detection properties of xenon gas, the analysis techniques 
described here have a direct application toward both low-level monitoring and background mea- 
surement tools for NEXT- 100. Concerning monitoring of key operational parameters, we have 
discussed how alpha particles can be used to track time variations at the percent level or better in 
the relative electroluminescent gain and, especially, in the electron drift velocity. Also, radon gas 
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diffusing in the detector is a potential background to the neutrinoless double beta decay search, see 
for example flZ7H . The work presented here is a first step toward establishing techniques to directly 
measure the concentration of radon gas in the detector. 
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